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Abstract 

In a chiral quark model where chiral symmetry is introduced via a non-linear a model, we 
evaluate the axial couplings 5u(0), (0) and 5^4 (0) related to the spin structure of the nucleon. 
Our calculation includes one-body and two-body axial current and pion absorption operators, 
which satisfy the Partial Conservation of Axial Current (PCAC) condition. While 5u(0) is 
dominated by the one-body axial current we find significant corrections due to two-body axial 
exchange currents in g%(0) and g^(0). Interestingly, the axial current associated with gluon 
exchange reduces g%(0) from 1 to 0.6. Our analysis shows that the so-called "proton spin crisis" 
can be resolved in a constituent quark model in which PCAC is satisfied. Furthermore, we use 
the PCAC constraint in order to determine the couplings of the r\ and r\' mesons to nucleons. 
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I. INTRODUCTION 



In a recent work [1] we have evaluated the axial form factors of the nucleon, g^, g n f n ~' po e 
and #t in the context of a chiral constituent quark model. Chiral symmetry was introduced 
via a linear a-model with pions and cr-meson degrees of freedom. The model is consistent with 
PCAC, up to some order in the nonrelativistic expansion of the axial quark operators. The 
PCAC constraint for the non-pole part of the axial current reads 

[H, A non - pole '°(q)] - q • A n "°' e (q) = % ^2 f n M 7r (q), (1) 

where H is the Hamiltonian of the three-quark system, including the center of mass motion, 
q is the three-momentum transfer of the weak gauge boson, A non ~ pole ^ (/i = 0, 1,2,3) stands 
for the non-pole part of the axial current operator, M n is the pion absorption operator, and 
f w = 92.4 MeV is the empirical pion decay constant. Eq.(l) plays a similar role in constraining 
the form of the axial currents as the continuity equation for the electromagnetic current. 

One of the main consequences of the PCAC relation for the axial current is that axial 
coupling of the constituent quarks, gA g , is not equal to unity. Instead, it is related to the 
pion-quark coupling constant, g nq , via a Goldberger-Treiman relation. The second important 
consequence of the PCAC constraint is the necessity to include axial exchange currents into the 
quark model. The presence of two-body potentials in H should be accompanied by two-body 
axial exchange currents and pion absorption operators if PCAC as formulated in Eq.(l) is to 
hold. The results of Ref. [1] show that the isovector axial coupling <?a(0) of the nucleon is 
dominated by the one-body axial current. Gluon, pion, and scalar axial exchange currents, 
although they are individually quite large, add to an overall correction of less than 2%. 

In the present paper we are interested in the two remaining axial nucleon couplings, namely 
the flavor octet isosinglet coupling g\(0) and its flavor singlet counterpart g^(0) related to 
the spin content of the nucleon. The spin fractions Aq carried by the quarks of flavor q are 
axial current matrix elements and can be evaluated in a constituent quark model. In the naive 
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nonrelativistic quark model, which uses only one-body axial currents, one obtains Au = 4/3, 
Ad = — 1/3 and As = and from them <?a(0) = 5/3 and g\(0) = <7^(0) = 1. On the other hand, 
the experimental numbers are: Au = 0.84 ±0.02, Ad = -0.42 ±0.02, and As = -0.09 ±0.02 [3] 
leading to g A (0) = 1.26 ± 0.04, g\(0) = 0.60 ± 0.08 and g\(0) = 0.33 ± 0.06. 

This marked disagreement between theory and experiment has often been interpreted as a 
severe shortcoming of the constituent quark model ("spin crisis"). Here, we show that once 
the PCAC violating and commonly used approximations: (i) axial charge of the constituent 
quarks g Aq is equal to unity i.e., the same as for QCD quarks, (ii) neglect of two-body axial 
exchange currents, are removed [1], the calculated spin fractions are in good agreement with 
the experimental numbers. 

The paper is organized as follows. In sect. II we review the relation between the axial 
current matrix elements and the spin fraction carried by the quarks as determined in deep 
inelastic lepton-nucleon scattering and octet baryon decays. Sect. Ill, reviews the non-linear 
(X-model of Manohar and Georgi, which is the basis of the chiral quark model used in sect. IV. 
The axial current operators are listed in sect. V, and our results for the nucleon spin structure 
are presented in sect. VI. There we also calculate the r/NN and rj'NN coupling constants. A 
summary of our results is presented in sect. VII. 

II. SPIN STRUCTURE OF THE NUCLEON 

In naive nonrelativistic constituent quark models, in which the three constituent quarks are 
in relative S-wave states, the nucleon spin is accounted for by the sum of the constituent quark 
spins. Since the late 1980's we know from deep inelastic scattering experiments that only a 
small fraction of the proton spin is carried by the spin of the QCD (or current) quarks, and 
that there is a nonzero contribution to the nucleon spin arising from strange quark- ant iquark 
pairs in the Dirac sea [2]. Recent experiments conclude that approximately 1/3 of the spin of 
the nucleon is carried by QCD quarks [3]. Calling Aq the fraction of the proton spin carried 
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by quarks of flavor q, we have from Ref. [3] 

Au = 0.84 ±0.02, Ad = -0.42 ±0.02, As = -0.09 ± 0.02. (2) 

In a parton model description of the nucleon, the spin fractions of the total nucleon spin 
carried by the individual quark flavors are given as 

Aq = [ dx(q + (x) + q + (x)—q-(x) — q~-(x)). (3) 
Jo 

Here, q±(x) is the probability of finding in the nucleon a QCD quark of flavor q with momentum 
fraction x of the total nucleon momentum having its helicity parallel (+) or antiparallel (-) to the 
proton helicity. The quark momentum distribution functions q±(x) can be further decomposed 

as 

q±(x) = ql al (x) + qr(x), (4) 

where q™ l (x) (q± a (x)) is the contribution of the QCD valence (sea) quarks; q±(x) stands for 
the antiquark distribution function in the sea. 

In deep inelastic lepton-proton scattering one is measuring the spin-dependent structure 
function g[{x, Q 2 ) of the proton at high momentum transfers Q 2 . The integral of this function 
rf = Jq g{(x, Q 2 ) dx is given in the scaling limit by the parton model result [4] 

T{ = -WA(z= — Au+— Ad+— As (5) 
1 2Y !8 18 18 

where e q is the quark charge of flavor q. The spin fractions can be expressed as axial vector 
current matrix elements 

(P,s\q-f^ 5 q\p,s) = Aqs^, (6) 

where is the spin vector. 

Additional information on the spin fractions can be extracted from the weak semileptonic 
decay of octet baryons. From neutron /3-decay one can extract the combination g\(0) = 
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Au — Ad which is nothing but the nucleon axial coupling constant <?a(0). Also, from SU(3) 
flavor symmetry and experimentally known hyperon /3-decays one obtains the spin fraction 
combination g\(0) = Au+Ad—2As. The two constants su(0) and <7^(0) govern all octet baryon 
/9-decays in the flavor symmetry limit. The combined DIS and hyperon /3-decay data give values 
for each Aq as quoted above. While the use of SU(3) symmetry is not without controversy and 
has been sometimes severely critizied [5], some authors think that SU(3) symmetry breaking 
effects are small, of the order of 10% [6]. 

If the Dirac sea were SU(3) flavor symmetric, g^(0) would have the interpretation as the 
fraction of the proton spin carried by valence QCD quarks, whereas the total fraction of the 
proton spin carried by all QCD quarks is given by the quantity 5^(0) = Au + Ad + As. Using 
the numbers [3] given above one obtains g° (0) = 0.33 ± 0.06, definitely smaller than 1, which 
is the expectation according to the naive nonrelativistic quark model. This is the origin of 
the so called "spin crisis". The experimental result may change somewhat. From the use of 
the flavor tagging technique at the HERMES experiment it would be possible to obtain the 
spin fractions directly from deep inelastic scattering experiments without relying on the SU(3) 
flavor symmetry assumption. Preliminary results from the experiment show that As > in 
disagreement with the previous determination [7]. This latter result has been critizied on the 
account that it would take <7^(0) out of any reasonable bounds [8]. 

Before we calculate the axial nucleon couplings we will review the SU(3) extended version 
of the non-linear a-model, which is appropriate for the evaluation of g\(0) and g° (0). 

III. THE NON-LINEAR a MODEL 

Let us start with the effective Lagrangian for the non-linear a model, which lends support 
to the notion of constituent quarks interacting via gluon and pseudoscalar boson exchange [9] . 
The effective Lagrangian is given by 

C = Hi'fip^ + - Mtttt + g Aq ^ l5 A^ + f -ftr j>E^E) -±tr (F^ F>») . (7) 
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Here, \1/ = (u, d, s) is the quark field and M = m q l is the quark mass matrix prior to any 
explicit SU(3) breaking. Furthermore, is the color SU(3)c covariant derivative 

D, = d, + igG„ G, = \ \) ■ G c ^, (8) 

where G M is the gluon field and the A^ are the color SU(3)c Gell-Mann matrices. F^ v is the 
gluon field strength tensor 

F, v = d,G v -d v G, + ig [G», G v \. (9) 

The fields V M and A^ appearing in Eq.(7) are given in terms of the matrix field £ derived from 
an octet of pseudoscalar boson fields <& 

v„ = \{id,e + £ f <m a, = l -{id,e - e%o- (io) 

The pseudo-scalar bosons <& are the Goldstone bosons (GB) associated with the spontaneous 
breaking of the SU (3)y x SU (3) a chiral symmetry. Following Ref. [9] we describe their dynamics 
in terms of a 3 x 3 matrix field S given by 

Y, = e iX *l F = e (11) 

where A is a vector whose components Xj are the eight Gell-Mann matrices of SU(3) flavor 
symmetry, and 4> stands for the eight GB fields <&_,-, j — 1 • • • 8. Here, i — 1,2,3 stands for the 
isotriplet n, % — 4, • • • 7 for the two doublets of K, and % = 8 for the rj. With this normalization 
F is the GB decay constant. Working in lowest order as we shall do 

F = f n = 92.4 MeV. (12) 

The matrix field 

^ = giA-S/2/. (13) 
transforms under S77(3)y x SU(3)a as 
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SU(3) V , SU(3) A , , , 

where V = exp(— iX ■ ay/2) and A = exp(— iX ■ a^/2) are global transformations belonging to 
SU(3)v and SU(3)a respectively. As for U, it is a unitary matrix field that depends on the 
axial transformation A and the $ Goldstone boson fields. 
On the other hand the quark field, transforms as 

SU(3) V SU(3)a 

$4>f — > A f/^r (15) 

while the gluon fields are invariant under both SU(3)y and SU(3)a transformations. 

At this point a comment on the value of g^ q in the effective Lagrangian is in order. Each 
term in Eq.(7) is separately invariant under the chiral group SU(3)yX SU(3)^. This is why one 
can introduce an axial quark coupling g^ q 7^ 1 in the model without violating chiral symmetry. 
The fact that gA q for constituent quarks can be different from unity is well founded, both from 
the phenomenological as well as the theoretical point of view. Weinberg [10,11] has shown that 
while constituent quarks have no anomalous magnetic moments, their axial coupling may be 
considerably renormalized by the strong interactions. Explicit calculation in [11] gives 

9M=l- g^j, (16) 

which corresponds to a 8% reduction of gAq- Using Witten's large Nc counting rules, one finds, 
that in contrast to the magnetic moment of the quarks, corrections to gA q appear already in 
order Nq [12]. Finally, investigation in the constituent quark structure in the Nambu and 
Jona-Lasinio model of Ref. [13] shows that gAq ~ 0.78. 

Explicit SU(3) breaking mass terms are introduced in the model as 

-vj7 e -ij 5 X-$/2f n Mq e -i 76 A-*/2/ 7r ^ ^ 

where M = diag(m^, m% m° s ) is the current quark mass matrix. After an expansion in powers 
of 1/ f n we get to dominant order 



M * + * 75 {A,-, M }^ -L * { Aj , {A fc , M }}M/ (18) 

The last term generates a mass term for the GB when the quark bilinear is replaced by its 
vacuum expectation value, see for instance Ref . [14] . 

When putting everything together with the Lagrangian in Eq.(7) we get to dominant order 

C = *(i 7 ^ — M — M )^ + l - dPQjdfij - l - m|.$5 - l -tr (F^ F^) 

/,-H *yA^ ^9^, + g^* {A* {A fe , M }}* *,.$ fc , (19) 

where M+M gives our constituent quark mass. The f jkl are the SU(3) structure constants and 
the are the masses of the GB. Expressions for these masses in terms of the current quark 
masses and the vacuum expectation value of the quark bilinear are easily obtained from Eq.(18) 
and are given in Ref. [14]. We assume that the mass terms break SU(3) flavor symmetry while 
preserving SU(2)j isospin symmetry. This means we shall take m" = m d . For this quantity one 
can use the value m|J = m d = m ~ 6 MeV, while m° can be determined now from the pion and 
kaon mass to be m° = ^f- 2m — m ~ 146 MeV. The actual values are not so important in this 
calculation, where the total constituent quark mass is the relevant quantity. 

In Eq.(19) we have apart from the kinetic energy and mass terms the gluon-quark interaction 
and GB-quark interaction terms. These interaction Lagrangians provide the vertices needed 
to build the one-gluon exchange potential of Fig. la, the one-boson exchange potentials of 
Fig. lb, the one-body contribution to the GB absorption operator of Fig. 3a, the two-body 
contribution to the GB absorption operator due to gluon exchange of Fig. 3b, and the two-body 
GB absorption operator due to one-boson exchange of Fig. 3d. In addition, the diagram in Fig. 
3c corresponding to two-body GB absorption due to the confinement interaction is obtained 
after the confinement potential of Fig. lc is introduced in sect. IV. 

In contrast to the linear a model, the coupling between the u and d quarks and the Goldstone 
bosons is now predominantly pseudo- vector in nature. Neglecting the pseudo-scalar term whose 
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coupling constant is given by m// x < 1 we obtain for the SU(2)j sector 

Ar* = ^l^T^d^, (20) 

with r being the isospin Pauli matrix and ir the pion field. Writing the coupling constant in 
the usual way as g nq /2m q we recover the Goldberger-Treiman relation at the quark level for 
the u and d quarks 

9M = fM, (21) 
m q 

where m q is the u and d constituent quark mass for which we use a value of m q = 313 MeV. The 
way g nq is fixed in Ref. [15] is unchanged here because for on-shell quarks the pseudo-vector and 
the pseudo-scalar quark-pion couplings are equivalent. This means we recover our result [1] 

g Aq = 0.774. 



A. Vector and axial currents 

The vector currents that derive from the Lagrangian in Eq.(19) and the transformation 
properties of the fields are given by 

Vf = + f jkl *y 75 A,* + f jkl Qkd^i (22) 

Using the equations of motion of the fields one finds that to dominant order the divergence of 
the vector currents are 

d^ = --n\ v M^. (23) 

This shows that only the vector currents corresponding to j = 1,2,3 (isospin) and j — 8 
(hypercharge) are conserved to that order. 

From these vector currents we can build up the electromagnetic current. With the quark 
charge given as 
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• (iA 3 + ^A 8 j (2 1) 



2 6 

the electromagnetic current will be 

J" = e ( V B+^vfj, (25) 

which is exactly conserved to all orders. 

Turning now to the axial currents one gets for the non-pole part of those currents 

A non- P ole,» = ^^^^ + _1 ^ ^ (2g) 

V2 V2 /tt 

The boson pole part of the axial current A^ B ~ pole '^ — y/^f^d^^j, which is completely deter- 
mined by the Goldstone boson fields, is omitted. 

The first term in Eq.(26) provides the extra vertex needed for the one-body (Fig. 2a), 
the two-body gluon (Fig. 2b), and two-body confinement (Fig. 2c) axial exchange current 
operators, and the second term in A non ~ pole ' M gives the one for the two-body axial current 
operator due to one-boson exchange (Fig. 2d). 

The divergence of the non-pole part of the axial current is given to dominant order by 

^ A non- P ole,u = _^ ^ ^ + ^ $ . + _J_ ^ ^ (2?) 

where once again the equation of motion has been used. The second term in Eq.(27) comes 
from the pseudoscalar quark-GB coupling that we shall neglect for the u and d sector. In its 
absence we have 

d»A n ° n - pole ^ = -V2U{d^ + m%)<$> 3 ■ j = 1,2,3,8. (28) 

This means we shall have exact PCAC for each term in the axial current, both at the relativistic 
level, and to the order we work in the expansion in f n . 

IV. THE NONRELATIVISTIC CHIRAL QUARK MODEL 



In this section we will briefly discuss the Hamiltonian of the chiral quark model, as motivated 
by the non-linear cr-model of the preceding section, and the wave function of the nucleon. 
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A. The Hamiltonian 



We are interested only in non-strange baryons so that only u and d quarks are needed. The 
Hamiltonian for the internal motion of the three quarks in a baryon is given by: 

h = E (™. + 4) " £ + § ((^ + <n* + ^ + ivm) >*) ■ (29) 

is the momentum operator of the j-th quark, and P is the momentum of the center of mass 
of the three quark system whose contribution to the kinetic energy is subtracted from the total 
Hamiltonian. Apart from the confinement potential (v con ^, the Hamiltonian includes two- 
body potentials from one-gluon (V 9 ), one-pion (V n ), and one- 773, (or $ 8 ), (V m ) exchange. The 
latter potentials are obtained from the Feynman diagrams in Fig.l. For the it and i] 8 meson 
exchange potentials we introduce a short distance regulator by means of the static vertex form 
factor 

/ A2 \ 1/2 

F < k H^n?) • < 3 °> 

Here, k is the three-momentum of the exchanged meson and A is the cut-off parameter. In 
coordinate space this leads to a very simple form for the potential where a second Yukawa term 
with a fictitious meson mass A appears. The actual expressions are given by 

t J S Al^f T ' ' Tt { °' ' V "> ' V '° (— " — J (31) 

((1a\ 2 1 A 2 1 /f- m V8 r p- A 18 r \ 

^ = (m S K^K, 3 ^ ■ ^ <- ■ ^ (V - V) < 32 > 

where rj, Uj and Tj are the position, spin and isospin operators of the j-th quark. The relative 
coordinate r is given by r = Yj — rj, and r = |r|. The factor 1/3 in the rjg potential comes from 
the normalization of the SU(3) flavor matrix A 8 . We will take the cut-off parameters equal 

A„ 8 = A^ = A. (33) 

The value of the cut-off A is fixed in Ref. [15], by fitting the size of the qq component of the 
pion, assumed to be of 0.4 fm as given in Ref. [16]. A value of A = 4.2 fm -1 results. For the 
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masses we use m n = 139 MeV, and assuming no octet-singlet mixing the Gell-Mann-Okubo 
mass formula gives m^ g = |(4m|- — m 2 .) from where m m ~ 565 MeV is obtained, which is very 
close to the mass of the physical 77 meson. 



The chiral quark model also contains the one-gluon exchange potential introduced by de 
Rujula et al. [17], and originally used to explain the spectrum of excited baryon states [18,19]. 
For this gluon exchange potential we use 




where a s is the effective quark-gluon coupling constant that we will take as a free parameter, A c 
are the SU(3) color matrices and v — r/r. Following the argumentation of Isgur and Karl [18] 
we do not consider the spin-orbit terms in the calculation because it is known that such terms 
would heavily distort the spectrum of the negative parity baryons. 



The constituent quarks are confined by a long-range, spin-independent, scalar two-body 
potential. For convenience a harmonic oscillator potential is often used 

(V™f). k = -a A< • K r 2 . (35) 

However, lattice calculations suggest that a linear radial function is more realistic. A linear 
confinement which a larger distances is screened by quark- ant iquark pair creation is found in 
some lattice calculations. It is known that a pure h.o. potential without any anharmonicity 
cannot reproduce the baryon mass spectrum [18,19], therefore we will use here a color screened 
confinement potential of the form 

(yconf^ ^ = _ a X c .\ c k (l- e'^) + C, (36) 

where a, \i and C are free parameters. 

B. The nucleon wave function 

Even though we work with the color screened confinement potential we will use the harmonic 
oscillator basis to expand our wave functions. We shall restrict ourselves to states up to N = 2 
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excitation quanta. The nucleon wave function is thus given by a superposition of five h.o. 
states: 

I AO = a Ss \S S ) + a s > s \S' S ) + a SM \S M ) + a Du \D M ) + a ?A \P A ) (37) 

A complete description of the wave functions used in Eq.(37) can be found in Ref. [19]. The 
mixing coefficients will be determined by diagonalization of the Hamiltonian in Eq.(29) in this 
restricted h.o. basis. 

The free parameters of the model are a s , a, /i, C and b, where b is the harmonic oscillator 
constant that appears in the radial part of the wave functions. The nucleon spectrum alone 
does not provide sufficient constraints in order to fix the parameters uniquely. For that reason 
we include also the A spectrum and the low energy electromagnetic nucleon properties, such as 
the magnetic moments and the charge and magnetic radii in order to fix the parameters. For 
details see Ref. [20]. The model parameters are shown in Table I while in Table II we give the 
admixture coefficients for the nucleon wave function. The Ss component clearly dominates. 

V. THE AXIAL CURRENT OPERATORS 

Using the vertices extracted from Eqs.(19, 26) it can be readily shown that each separate 
contribution in Figs. 2-3 satisfies in momentum space 

^ A non- P ole,» = • ^ f n M f B (38) 

for the components j = 1,2,3,8 of the current and for non-strange quarks. Here, q is the 
four-momentum transfer, y[™ on -P ole ' v j s ^he axial current amplitude corresponding to Fig.2(a- 
d), and Mf B the Goldstone boson absorption amplitude of Fig.3(a-d). This result is a direct 
consequence of Eq.(28). 

Let us now give the explicit expressions for the different axial current operators that cor- 
respond to the nonrelativistic reduction of the Feynman diagrams in Fig.2(a-d). Even though 
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we will only need the spatial part of these operators we give here, for completeness, also the 
expressions for the zeroth component, i.e., for the axial charge operators. In the rest of this sec- 
tion we drop the non-pole affix of the axial current operators in the understanding that we will 
always refer to its non-pole part. The one-body isovector axial current operators (j = 1,2,3) 
are 

r!?> ..... i 



k 

Aj, imp = g Aq ]T - *= a k . (39) 
k V2 



In order to generalize the isovector axial current to the flavor octet isosinglet axial current 
we replace in the axial vertex by \/3 A| to get for the j = 8 case 

<, imP = a Aq e eiq ' rfc ^r ° k ■ (q + 2pfc) 

A; 9 

A 8 , imp = E eiqrfc ( 40 ) 

where we have taken into account that nucleons are made of u and d valence quarks so that we 
only need the left upper corner of the SU(3) flavor matrix A 8 = l/V3diag(l : 1, —2) which is 
proportional to the unit matrix. 

As in the derivation of the two-body potentials, we consistently neglect all non-local terms 
in the exchange currents. For the one-gluon exchange operators we then obtain for j = 1, 2, 3 



3,9 ~ q hi^ m \ 1 IV2 



-i((7 fe • r) q 

+ (3(((7 fc + a,)-r) f - (C7 fc + (7,)) ^ 

+ (k^l)\ (41) 



+ y((7 fc +(7,) (J(r) 



and for the 8-th component 
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< 9 - ~9mE £^ • A? | (<r fe x a,) • r + (*<-► Z)}^ 



fc<Z 16 ™S 



-i(<r fc • r) q 



(3((<7 fe +a t ) ■ r) f - (o- fc + ^ ^ 
+ (fc<->Z)l. 



+ y(<7fc+<7,) (J(r) 



Now the confinement exchange contribution reads for j = 1, 2, 3 



A - n 



-%{a k • q) q + (q • V r ) a k - (a k ■ V r ) q 



+(£;^Z)j ^ con/ (r) 



and similarly for the 8-th component 



4° = n 

^conf — U 



A 8 , con/ EE^ A(? — e iq ' 

^'"g fc<Z I 



■i(a k • q) q + (q • V r ) 0fc - (a k • V r ) q 

+(*;<-> Z)j \/ con/ (r) . 



Finally axial pion-pair exchange current operators are given for j = 1, 2, 3 by 
„o _ ^ 1 A 2 



2/2 4vr A 2 - m 2 " g 

9Ag 1 1 A 2 



Tfc X J l)b) - ^i(r) a, • r + (Jfe <- I) 

V2 r 



A • 



2/2 2m g 4tt A 2 - m 2 



E 



-(Tfc X T;)k"> 
^2 



(o-fc x a t ) - ^i(r) 
r 



1 



+i(o-fc x q) (<7, • r)- ^i(r) 
- (07 • r) (<7 fe x r) \ y 2 {r) 



yi (r) and ( r ) are given by 
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y ± ( r ) = ml Yxim^r) - A 2 Yi(Ar) 
y 2 (r) = m 3 Y 2 (m n r) - A 3 F 2 (Ar) 



where 



= — (i + -) 

p-x o o 

^) = —(! + - + 9 



For the 8-th component we have for the use with non-strange baryons 

4U = (46) 
A 8)7r = 0. (47) 



VI. SPIN CONTENT OF THE NUCLEON 

Because the spin operator of Dirac particles is related to the axial current operator we 
can obtain information concerning the spin content of the nucleon from the axial form factors 
calculated in the present chiral quark model. As explained in sect. II, the spin fractions can be 
expressed in terms of matrix elements of the axial current operators. The latter are based on the 
larger SU(3)yX SU(3) J 4 symmetry underlying the Lagrangian of the non-linear a model, which 
contains SU(3) flavor symmetry as a subgroup. There are now three axial vector couplings of the 
nucleon: the flavor octet isovector axial coupling 5u(0), the flavor octet isosinglet (hypercharge) 
axial vector coupling <?^(0), and the flavor singlet axial coupling g° A (0). 

A. The flavor octet isovector axial coupling gA(0) 

The nucleon isovector axial coupling 5u(0) can be obtained from the spin quantities defined 
above as 

g A {0) = Au- Ad= 1.26 ±0.04. (48) 
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This value agrees within errors with the one measured in neutron /3-decay given by 
fiU^ P (0) = 1.2670(35). In Table III we give the results that we obtain in our model. The 
impulse contribution clearly dominates. Exchange currents account in this case for less than 
9% of the total value which is in reasonable agreement with the experimental number. As 
explained in Ref. [1] the fact that the axial coupling of the constituents quarks g Aq is smaller 
than unity is mainly responsible for the improvement of the calculated g A (0) compared to the 
naive quark model result of 5u(0) = 5/3. 



The flavor octet isospin singlet axial coupling g A (0) gives the fraction of the proton spin that 
is carried by valence QCD quarks if the Dirac sea were SU(3) flavor symmetric. Experimentally 
it is determined as 



We evaluate g A (0) using the 8-th component of our axial operators given in sect. V and 
the nucleon wave function given in sect. IV. From there we see that only impulse and gluon 
exchange currents contribute to g A (q 2 ) at four- momentum transfer q 2 = 0, which for nucleons 
also implies q 2 = 0. 

In the first line of table 4 we show our results with g Aq = 1- The impulse contribution is 
then very close to 1. This has to be so since A 8 ,i mp (q 2 = 0) is gA q times the total spin operator 
J2k a k and the quarks are mainly in relative S-wave states. 

Two-body exchange currents play a very important role for this observable. The axial 
gluon exchange current, which effectively describes the qq pairs in the nucleon reduces the 
value obtained in impulse approximation by some 40%. In addition, there is a reduction of 
g A (0) because g Aq < 1- The total result is shown in the second line in Table IV. Our total 
value is some 30% smaller than the present experimental data. However, because SU(3) flavor 



B. The flavor octet isospin singlet axial coupling g^ 



(0) 



g 8 A (0) = Au + Ad- 2As = 0.60 ± 0.08. 



(49) 
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symmetry is broken the condition g\ = g Aq need not be exact. Explicit calculation of axial 
constituent quark couplings in the Nambu-Jona Lasinio model shows that g\ q ^ gAq [21]. In 
summary, the experimentally observed spin fraction carried by the valence QCD quarks can be 
obtained if axial gluon exchange currents are included in the nonrelativistic quark model. 

C. The flavor singlet axial coupling g^(0) 

The other independent quantity that one can construct from the Aq's is g°(0) which is 
experimentally determined as 

g° A (0) = Au + Ad + As = 0.33 ± 0.06. (50) 

The singlet axial coupling <?^(0) is the nucleon spin fraction carried by all quarks, i.e., valence 
plus sea quarks. 

In order to evaluate fi^(0) in our model we have to extend chiral symmetry to the group 
U(3)v x U(3)a- The QCD Lagrangian is trivially invariant under U(l)v which is the symmetry 
associated with baryon number conservation. As for U(1)a even though the QCD Lagrangian 
is invariant in the chiral limit we know that this symmetry is broken by quantum effects. 
This is called the U(1) A anomaly [22]. If the anomaly were not present we would have another 
Goldstone boson $ i n the chiral symmetry limit m q — > 0. However, the presence of the anomaly 
makes m$ nonzero even in the chiral limit of zero current quark masses. The natural candidate 
for this extra boson is the physical rj' meson whose mass is much larger than the masses of the 
pseudoscalar octet mesons, which are the quasi-Goldstone bosons of chiral 577(3). In the large 
N c limit of QCD the U(1) A anomaly is absent. In that limit the massless QCD Lagrangian 
has the larger symmetry U(3) v x U(3) A and there are nine Goldstone bosons associated with 
the spontaneous chiral symmetry breaking to the subgroup U(3)v- 

In that large N c limit the Goldstone boson dynamics can be described in terms of the 
matrix field 
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where A = ^2/3/ and £ is given in Eq.(ll). The new boson field, $o, is invariant under 
6 r (3)y x S77(3)a whereas under U(1)a it transforms as 

$ ^ $ o + ao/ w - (52) 

Here, ao is the parameter associated with the U(1)a transformation. As for the other boson 
and quark fields the matrix field £ defined in Eq.(13) is invariant under U(l)yX U(l)^ and the 
quark field * transforms under U(l)y as 

being invariant under U(1)a because the U matrix field does not depend on a or $ . 

With these transformation properties, the most general Lagrangian invariant under £7(3) y x 
U(3)a has as the first few terms 

+|^ ^ 7 M 75 A ^ + \ d,^ d^ + t g°k^ d^o tr (Ec^t) (54) 

where g Aq and g° are free parameters not fixed by chiral symmetry requirements alone. 

Including also the current quark mass breaking term, where now the expression in Eq.(17) 
is changed to include a factor e ~ i75Ao * 0//2 ^' r to the left and to the right of M , and making an 
expansion in powers of 1/ f n we get to dominant order 

£ = ^{i^d^ - M - M )^ 

+ \d^ 3 d^ 3 - \ m% $J + \d»%d,% - \ m| o <&g - ml s o $ 8 $ 
-^H-f-ftX^dfij + f^^Ts XoVdpQo + ^7 5 {A Q , M }^ $ a 

^Jtt ^Jtt ^JlT 

- fjk&'fXj* *kd,*$i + HK, {A/?, M }} (55) 
where the a and (5 indices run from to 8. 
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The mass term for the <3> meson has to be supplemented by a quantity that is not vanishing 
in the chiral limit of zero quark masses, see for instance Ref. [23]. Due to the $ 8 and <3> mass 
mixing term those bosons are no longer the physical states. The physical states are the r\ and 
i]' mesons obtained by diagonalizing the mass matrix and that are given in terms of $8 an d $o 

as 

rj = cos#p $ 8 — sin#p $! 

rj = sin9p $ 8 + cosO P $i (56) 

where Op is the pseudoscalar mixing angle. There are different determinations of Op in 
the literature. In the theoretical analysis of Ref. [24] a value of 9 P = —12.3° is obtained 
while a phenomenological analysis of different experimental data done by the same authors 
gives 9p = —15.4 ± 1°. In Ref. [25] the study of vector-pseudoscalar-photon decays gives 
dp = — 17 ±2.9°. A recent evaluation in Lattice QCD done in Ref. [26] gives the result 
6 P = —10 ± 2°. Finally the value favored by the Particle Data Group [27] is P —20° which 
is obtained from the analysis of pseudo-scalars two-photon decays. Different determinations 
range from —10° to —20°. 

Turning now to currents it is easy to see that for the vector currents with components 
j = 1,2,3,8 there is no difference from the results stated in Eqs. (22,23). The new U{l)v 
symmetry has an associated vector current that is given by 

Kf = *Vy*- (57) 

Apart from a constant factor, this is the baryon number current and is exactly conserved 

0^ = 0. (58) 

For the axial currents, and compared to the results in sect. Ill, we have that the divergence 
of the j — 8 component is now given by 

dp A n 8 on - pole >» = -V2f n {cos Op (d, + mf} rj + sin P (d, + mj) r]'} (59) 
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while for the j = 1,2,3 components there is no change. As for the new axial current corre- 
sponding to the U(1)a symmetry we have for the non-pole part 

AT^* = gX.v-fTtfyf (60) 

whose divergence is given by 

d„ A n on ~ P ° le '>* = -V2f v {- sin 6 P &> + mj) 77 + cos 6 P (d, & + mj) 7/} ■ (61) 

The two PCAC relations in Eqs. (59,61) translate into the following relations in momentum 

space 



A non-pole, M = { ^ ^ j cos ^ p M V + sin 6 P M V ' J 



where j\ non -p° le ^ stands here for the non-pole part of the axial current amplitude and M n > v ' 
is the absorption amplitude for the 77,77' mesons. Relations in Eq.(62) are exactly satisfied for 
each independent contribution namely: impulse of Figs. 2-3a, one-gluon exchange of Figs. 2-3b 
and confinement exchange of Figs. 2-3c. 

To evaluate g° A (0) we will need the non-relativistic operators associated with Aq OTI po e ' M 
current. From Eq.(60) and our discussion in sect. Ill it is clear that there will be impulse, 
gluon and confinement exchange current contributions similar to the ones introduced in sect. V 
for the 8-th component of the axial current. In fact, with the correct normalization and in the 
case of non-strange baryons, the operators are exactly the same as those given in sect. V for 
the 8-th component of the axial current with the sole replacement of g Aq by g Aq . 

Our results for g A (0)/g Aq are the same as in Table IV for g A (0)/g Aq . In order to reproduce 
the experimental value of g A (0) = 0.33 we need 

9% = 0-553 (63) 
which is 30% smaller than g Aq . The ratio 
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C = QaJqa, = 0.714 (64) 

obtained here is at variance with other empirical determinations like \(\ = 1.16 in Ref. [28] 
or ( = —1.2 in Ref. [29]. It is however, in agreement with the result g Aq = 0.54 obtained in 
Ref. [21]. 

A criticism is in order here. In our Hamiltonian we have used an i] 8 , or $ 8 , meson exchange 
potential with a r/ 8 coupling to u and d quarks given by 5 , A g /(2/ 7r V / 3). With i] and r/ 8 mixing 
one would have to change the rjs exchange potential into corresponding r\ and rj' exchange 
potentials with couplings to u and d quarks given by 

tt^~ — 7= ( c °s Op QAa — sin dp y/2 q° An ^ 

2m q 2f w V3\ q ' 

l^r^ MF9M+coa0p ^ s% )- (65) 

While this change will affect the determination of parameters and the nucleon wave function, 
the fact that the 7]$ and the rj and meson masses are much larger than the pion mass makes 
the corresponding potentials not so relevant in this respect. Their effect can be compensated 
by very small changes in the free parameters of the model. 

D. Effective couplings of the i] and rj' mesons to nucleons 

Using PCAC at the baryon level one can get relations between <?^(0) and g A (0) on one side 
and the rj and rj' coupling to nucleons on the other. Those relations are 

9 v 'nn(0) 1 



9 A (0) = V3fA cos^^^ + sin^ 



M 



N 



(66) 



and similarly 



(o) 



^°>=^M- sin9p ^ +cosS ^Sn (67) 

from where the r]NN and rj'NN coupling constants can be obtained as 
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g vNN = {cos^^(O) - sinflp V2g° A (0)} 

QrfNN = -^j {sin0 P <&(0) + cos^p V2g° A (0)} . (68) 

Using our theoretical results for g%(0) and g A (0) (the second one is in fact fixed to the experi- 
mental value), and with the use of 9 P in the range [—20°, —10°] we get 

g V NN = 3.18 - 3.52 

g v , NN = 2.22 - 1.63 (69) 

or 

a ,, NN = = o.80 - 0.99 

47T 

a 2 

a v , NN = = 0.39 - 0.21. (70) 

4:71 

There are various determinations of o^jvjv an d ^'nn in the literature. Typical values for 
a V NN suggested by the Bonn model [30] are in the range 3 — 5, much larger that our evaluation. 
From a practical viewpoint one need not worry too much about this discrepancy as the rj meson 
is not very relevant for the fit of the nucleon-nucleon phase shifts. Also its contribution to the 
nuclear binding is very small for normal nuclear densities. Besides, in iViV potentials the r]NN 
coupling is usually treated as an effective parameter in which contributions due to two-meson 
exchanges are absorbed. The analysis of the — > r/n reaction in Ref. [31] results in values 
for a^NN in the range 0.6 — 1.7 in good agreement with our determination. Smaller values are 
found in Ref. [32] where the analysis of the ao — > r\ix decay gives ol^nn = 0.315. Also the use 
of QCD sum rules on the light cone [33] gives a v ^N = 0.35 ± 0.1. 

As for a v 'NN, a recent theoretical analysis of pp scattering near the production threshold [34] 
gives ay at at < 0.5. Another analysis [35], this time of the near threshold rj' photo-production, 
obtains ay at at ~ 0.22. Both results are in good agreement with our evaluation. 
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VII. SUMMARY 



We have investigated the axial form factors related to the spin structure of the nucleon in 
the framework of a chiral quark model where chiral symmetry is introduced via a non-linear 
a model. We have included one-body and two-body axial currents as required by the PCAC 
condition. For <yu(0) we obtain a value of 1.34 in good agreement with experiment. We find the 
5u(0) is dominated by the one-body axial current. For <7^(0), related to the spin fraction carried 
by the valence quarks (if the sea is flavor symmetric), we obtain 0.47. This surprising result 
within the context of a nonrelativistic quark potential model, where one would expect a value 
of 1, is due to two factors. First, the inclusion of the relativistic correction provided by the 
axial one-gluon exchange current, and second the fact that the constituent quark axial coupling 
constant gA q satisfies the Goldberger-Treiman relation of Eq.(21) and is therefore smaller than 
unity. Due to the internal structure of a constituent quark, also its coupling to the singlet axial 
current g° Aq gets renormalized. Its value is considerably smaller than unity. 

Using our predictions for g A (0) and g\ (0) as well as the PCAC constraint, we calculate the 
effective coupling constants of the rj and rj' mesons to nucleons. The values obtained for these 
coupling constants are in agreement with different theoretical analyses done by other groups. 
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TABLE I. Quark model parameters: b is the harmonic oscillator constant, a s is the quark-gluon 
coupling strength, a is the confinement strength, jj, the color screening length, and C a constant term 
in the color screening confinement potential. 



b[fm] 




a[MeV] v[fm- 1 ] C[MeV] 


0.591 


0.987 


244.041 0.770 -279.881 


TABLE II. 
of Eq.(29) in 


Admixture coefficients in the 
an N = 2 harmonic oscillator 


wave function of Eq.(37), evaluated with the Hamiltonian 
space. 


a S s 


a S ' s 


a S M a D M a P A 


0.9585 


-0.1475 


-0.2344 -0.0672 0.0011 



TABLE III. Results for #a(0). Imp., impulse contribution ; g, gluon exchange contribution ; tt, 
pion exchange contribution ; Conf., confinement contribution; Total, total result. 

Imp. g 7r Conf. Total Exp. 



g A (0) 1.228 -0.206 0.321 1.343 1.2670(35) 



27 



TABLE IV. Results for g\(0) and g A (0). Imp., impulse contribution; g, gluon exchange contribu- 
tion; Total, total result. The constituent quark axial couplings are given by g Aq = g\ q = 0.774 and 
9 % = 0-553. 

Imp. g Total Exp. 
g%(0)/g Aq 0.9909 -0.3861 0.6048 

g\(0) 0.7670 -0.2988 0.4682 0.60(8) 
g° A (0)/g° Aq 0.9909 -0.3861 0.6048 

g° A (0) 0.5480 -0.2135 0.3345 0.33(6) 
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FIG. 1. Feynman diagrams for two-body exchange potentials, a) One-gluon exchange potential, b) 
One-pion and one-r/g exchange potential, c) Confinement potential. 
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FIG. 2. Feynman diagrams for the axial operators, a) Impulse, b) One-gluon exchange, c) Con- 
finement interaction, d) Pion exchange. 
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FIG. 3. Feynman diagrams for the Goldstone boson (GB) absorption operators, a) Impulse, b) 
One-gluon exchange, c) Confinement exchange, d) Pion exchange. 
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